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It has often been suggested that the cosmic positron excess observed by the HEAT experiment could 
be the consequence of supersymmetric dark matter annihilating in the galactic halo. Although it is 
well known that evenly distributed dark matter cannot account for the observed excess, if substantial 
amounts of local dark matter substructure are present, the positron flux would be enhanced, perhaps 
to the observed magnitude. In this paper, we attempt to identify the nature of the substructure 
required to match the HEAT data, including the location, size and density of any local dark matter 
clump(s). Additionally, we attempt to assess the probability of such substructure being present. 
We find that if the current density of neutralino dark matter is the result of thermal production, 
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10 or less) conditions must be present in local substructure to account for the 
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I. INTRODUCTION 

A large body of evidence has accumulated in the fa- 
vor of cold dark matter. This body of evidence includes 
observations of galactic clusters and large scale structure 
[1] , supernovae [2] and the cosmic microwave background 
(CMB) anisotropies [3,4]. At the 2er confidence level, the 
density of non-baryonic [5], and cold, dark matter is now 
known to be VL CGM h 2 = 0.113 j^gjj [4]. 

A compelling dark matter candidate is provided by su- 
persymmetry [6]. In supersymmetric models which con- 
serve R-parity [7], the lightest supersymmetric particle 
(LSP), is stable. Furthermore, in many supersymmetric 
models, the LSP is the lightest neutralino, a mixture of 
the superpartners of the photon, Z and neutral Higgs 
bosons, and is electrically neutral, colorless and, there- 
fore, a viable dark matter candidate. If such a parti- 
cle were in equilibrium with photons in the early uni- 
verse, as the temperature decreased, a freeze-out would 
occur leaving a thermal relic density. The temperature 
at which this occurs, and the density which remains, de- 
pends on the annihilation cross section and mass of the 
lightest neutralino. Supersymmetry is capable of provid- 
ing a dark matter candidate with a present abundance 
consistent with WMAP and other experiments. 

Many methods have been proposed to search for ev- 
idence of supersymmetric dark matter. These include 
experiments which hope to measure the recoil of dark 
matter particles elastically scattering off of a detector (di- 
rect searches) [8] , experiments which hope to observe the 
products of dark matter annihilation (indirect searches) 
and, of course, collider experiments [9]. Indirect searches 
include searches for neutrinos [10], gamma-rays [11], anti- 
protons [12] and positrons [13,15]. 

In 1994 and 1995, the High-Energy Antimatter Tele- 
scope (HEAT) observed a flux of cosmic positrons well in 
excess of the predicted rate, peaking around ~ 10 GeV 
[16]. This result was confirmed by another HEAT flight 



in 2000 [17,18]. Although the source of these positrons 
is not known, it has been suggested in numerous pub- 
lications that this signal could be the product of dark 
matter annihilations, particularly within the context of 
supersymmetry [13,15]. 

If the dark matter is evenly distributed in our local 
region (within a few kpc), the rate of annihilations will 
be insufficient to produce the observed excess. It has 
been suggested, however, that if sufficient clumping were 
present in the galactic halo, the rate at which such par- 
ticles annihilate could be enhanced enough to accommo- 
date the data. At this time, we have very little infor- 
mation regarding the presence of any dark matter sub- 
structure in our local region. This makes any predictions 
of the positron flux very difficult to make. With sophis- 
ticated numerical simulations and analytical models of 
galactic substructure, however, such predictions can be 
made on a statistical basis. In this article, we attempt 
to discern whether it is possible for supersymmetry to 
provide a flux of positrons sufficient to account for the 
observed excess, and if so, how much local dark matter 
substructure would be needed and what is the probability 
that the required amount of local substructure is present. 
We find that, although such an explanation is possible, 
the probability of the necessary local substructure being 
present is very small. 



II. POSITRONS FROM DARK MATTER 
ANNIHILATIONS 

Positrons can be produced in several neutralino an- 
nihilation modes. For example, they can result from 
the decay of gauge bosons produced in the interactions 
X°X° ~ ¥ ZZ or x°X° W + W~ , producing positrons 
of energy ~ m x o/2. A continuum of positrons, extend- 
ing to much lower energies, can also be produced in the 
cascades of particles produced in annihilations includ- 
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ing fermions, Higgs bosons and gauge bosons. The spec- 
trum of positrons produced in neutralino annihilations 
can vary significantly depending on the mass and anni- 
hilation modes of the LSP. 

If the lightest neutralino is lighter than the W and Z 
bosons, their annihilations will be dominated by x°X° — * 
bb with a small (typically a few percent) contribution of 
X°X° ~ * t + t~ . Assuming that these annihilations are 
dominated by bottom quarks, the spectrum of positrons 
produced will depend only on the mass of the LSP. If the 
LSP is heavier, the annihilation products can be more 
complicated, often dominated by several modes includ- 
ing xx ~ * W + W~, x°X° ~ * ZZ or X°X° —> tt as well as 
X°X° ~> bb and x°X° ~ * t + t~. Also, a small fraction of 
annihilations can directly produce an electron-positron 
pair, x°X° ~ * e + e~, although this occurs very rarely, 
making its impact generally negligible. In our calcula- 
tions, the spectrum from each annihilation mode, includ- 
ing cascading, is calculated using PYTHIA [19], as it is 
implemented in the DarkSusy package [20]. 

In figure 1, we show the positron spectrum from neu- 
tralino annihilations for the most important annihilation 
modes. Solid lines represent the positron spectrum, per 
annihilation, for x°X° bb, for LSPs with masses of 50, 
150 and 600 GeV. The dotted lines are the same, but 
from the process x°X° ~ > t + t~ . Gaugino-like annihila- 
tions typically produce a spectrum which is dominated 
by bb at low energies, with contributions from t + t~ only 
becoming important at energies of about half the LSP 
mass and above. 
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FIG. 1. The positron spectrum from neutralino annihi- 
lations for the most important annihilation modes. Solid 
lines represent the positron spectrum, per annihilation, for 
X°X° -> bb, for LSPs with masses of 50, 150 and 600 
GeV. The dotted lines are the same, but from the process 
X°x" —> t + t~ . Dashed lines represent positrons from the pro- 
cess XX — * W + W~ for LSPs with masses of 150 and 600 GeV. 
The spectrum from XX ~ * ZZ is very similar. The positron 
spectrum from WIMP annihilations will be the weighted sum 
of the spectra over the annihilation modes, such as those 
shown here. See text for more details. 



For neutralinos with a substantial higgsino component, 
annihilations to gauge boson will often dominate. Dashed 
lines represent positrons from the process XX ~~ > W + W~ 
for LSPs with masses of 150 and 600 GeV. The spectrum 
from XX ~ * ZZ is very similar. Notice that the positron 
spectrum via gauge bosons is slightly harder than from 
bb. 



III. POSITRON PROPAGATION 

We use a standard diffusion model to calculate the ef- 
fect of propagation on the observed positron spectrum. 
In this model, positrons, being charged particles, move 
under the influence of interstellar magnetic fields. The 
galactic magnetic field is tangled, resulting in motion 
which is well described as a random walk. During this 
motion, positrons lose energy via inverse Compton and 
synchrotron processes. 

The diffusion-loss equation describing this process is 
given by 
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where dn e + / dE e + is the number density of positrons per 
unit energy, K(E e + , x) is the diffusion constant, b(E e + , x) 
is the rate of energy loss and Q(E e + , x) is the source term. 

We parameterize the diffusion constant [21] and rate 
of energy loss by 



K(E e +) = 3 x 10 
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+ E°J 



cm 2 s 



and 



b(E e +) = lir w E 2 e 



(2) 



(3) 



respectively. b(E e +) is the result of inverse Compton 
scattering on both starlight and the cosmic microwave 
background [22]. The diffusion parameters are con- 
strained from analyzing stable nuclei in cosmic rays (pri- 
marily by fitting the boron to carbon ratio) [23]. 

In equations 2 and 3, we have dropped the dependence 
on location. We treat these as constant within the diffu- 
sion zone. For the diffusion zone, we consider a slab of 
thickness 2L, where L is chosen to be 4 kpc, the best fit 
to observations [21,23]. The radius of the slab is unim- 
portant, as it is larger than the distances which positrons 
can propagate at these energies. Outside of the diffusion 
zone, we drop the positron density to zero (free escape 
boundary conditions). 

We solve equation 1 following the procedure described 
in Ref. [15]. For detailed descriptions of two zone diffu- 
sion models, see Ref. [15,23,24]. 
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IV. SUBSTRUCTURE 

Given the shape of the primordial power spectrum for 
CDM fluctuations, dark matter halos are expected to 
form hierarchically from the merging of smaller bound 
systems, the cores of which survive as dense substruc- 
ture in present-day halos. This merger process has been 
studied extensively using numerical TV-body simulations. 
While simulations should provide an accurate picture of 
halo formation through merging, they cannot easily track 
the evolution of substructure in the densest regions of ha- 
los. Even the highest resolution simulations show little 
or no substructure in their innermost regions, but this 
is probably due to numerical effects rather than physical 
ones. To avoid underestimating the number of dark mat- 
ter clumps that might survive in the solar neighbourhood, 
we will use a recently developed semi-analytic model of 
halo formation [32]. The model treats the evolution of 
dark matter lumps in a way that is less sensitive to their 
mass or the mean density of their environment. Com- 
parisons with simulations suggest that if anything, this 
model may overestimate the amount of substructure in 
halos [32]. 

We will outline the basic features of this model here; 
for detailed information, see Ref. [32]. The model gener- 
ates sets of random but representative merger histories 
for dark matter halos, using the merger-tree algorithm of 
Ref. [28]. The merger trees are 'pruned' [32] to determine 
whether halos merging with the main system contribute a 
single lump or several lumps (corresponding to their own 
undigested substructure) to the substructure of the main 
system. The dynamical evolution of individual lumps 
is determined using the analytic model of satellite dy- 
namics developed in Ref. [29], which includes the effects 
of dynamical friction, tidal mass loss and tidal shock- 
heating. The background mass distribution is taken to 
be the sum of a spherical dark matter halo (with a cos- 
mological density profile), a spherical bulge component 
and an exponential disk. The dark matter halo grows in 
mass according to its merger history, and changes in con- 
centration following the relations proposed in Ref. [25]. 
The growth of the baryonic components is modelled as 
described in Ref. [27] in order to produce present-day 
systems with roughly the properties of the Milky Way. 

We run this model using a large set of input merger 
trees with a mass resolution of 5 x 10 7 Af Q , as well as a 
smaller number of trees with a resolution of 1 x 10 7 Mq , in 
order to test for convergence in our results. Within each 
model halo, we then select every surviving lump in the 
solar neighbourhood (7 kpc < R < 10 kpc). Given the 
average amount of stripping experienced by these sys- 
tems, our results for this region are complete down to 
roughly 3 x 1O 6 M and 6 x 1O 5 M for the low-resolution 
and high-resolution models, respectively. For each of 
these systems, we calculate a quantity proportional to 



the annihilation rate, / x M^ lump /T4i ump , assuming they 
start out with an NFW density profile with a concentra- 
tion given by the ENS concentration relations [25], and 
that stripping changes their density profile as described 
in Ref. [26], and the value of / as described in [30]. Here, 
/ is defined as J(p 2 /p 2 )dV. 

We note that assuming a density profile with a central 
slope close to r~ 15 for the subhalos would result in anni- 
hilation rates up to 20 times larger [30] , but a profile this 
steep is disfavoured by the most recent numerical results 
[31]. 
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FIG. 2. The number density of dark matter clumps be- 
tween 7 and 10 kpc from the galactic center, as a func- 
tion of the minimum value of / x Af^ lump /V c i ump considered. 
The thick solid line shows the results for merger trees com- 
plete down to ~ 3 x lO 6 Af0, while the thick dashed line 
shows the results for higher-resolution trees complete down 
to ~ 6 x 1O 5 M0. The thick dotted line shows results for a 
different disruption efficiency (see text). The thin lines show 
the ±l-<7 halo-to-halo variation. 

Figure 2 shows the number density of dark matter 
clumps between 7 and 10 kpc from the galactic center, 
as a function of the minimum value of / x -M^iump/^ciump 
considered. The solid lines show the results for the 
low-resolution trees (the thick line indicates the average, 
and the thin lines show the 1-er halo-to- halo variation), 
while the dashed lines show the results for the high- 
resolution trees. The two agree to within a fraction of 
the halo-to-halo scatter, for values of / x M 2 lnmp /V c i ump 
above their respective resolutions, indicating our esti- 
mate is reasonably independent of merger-tree resolution. 
Clearly the halo-to-halo scatter is very large; some sys- 
tems have almost no clumps close to the solar radius, for 
instance. There is a fairly well-defined upper envelope 
to the distribution, however. For example, for values of 
/ x ^lump/^ciump > 2 x 10 7 M|pc~ 3 , no more than 
about 10~ 3 clumps per cubic kpc are predicted. Further- 
more, we have reason to believe this is a very conservative 
upper limit. Not only do simulations find less substruc- 
ture locally, as mentioned previously, but even within 
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the semi-analytic model there is some uncertainty as to 
whether this many systems really survive. We have as- 
sumed that subhalos on general orbits can survive strip- 
ping down to one tenth of the critical radius defined by 
Ref. [26] for circular orbits, by which point they have typ- 
ically lost more than 99% of their original mass ('model 
B' in Ref. [32]). If we make the less extreme assumption 
that systems are disrupted once they have been stripped 
down to half their critical radius (97-98% mass loss), 
then far fewer systems survive in the solar neighbour- 
hood, as indicated by the dotted line in figure 2. 



V. POSITRON SPECTRA FROM DARK 
SUBSTRUCTURE 

We will now consider a single clump of dark matter 
particles. We justify this assumption by pointing out 
that the probability of having a single clump within a 
few kpc which contributes some particular amount to the 
annihilation rate is considerably larger than the proba- 
bility of having multiple clumps which sum to the same 
contribution. This can be observed in figure 2. 

The annihilation rate in a dark matter clump is de- 
pendent on the clump's mass and density profile as well 
as the dark matter particle's mass and annihilation cross 
section. Treating the effect of these quantities as a free 
parameter used to normalize the flux of positrons to ob- 
servations, we can assess if the observed spectrum can be 
fit by such a scenario. 

We will consider five representative supersymmetry 
scenarios. First, a light (50 GeV) neutralino which anni- 
hilates 96% of the time into bb and 4% to t + t~. Such a 
particle could be gaugino-like or higgsino-like, as below 
the gauge boson masses, these modes dominate for either 
case. 

Second, we consider two cases for an intermediate mass 
(150 GeV) neutralino. One which annihilates as de- 
scribed in the previous case, and another which anni- 
hilates to gauge bosons (58% to W + W~ and 42% to 
ZZ). These neutralinos are typically gaugino-like and 
higgsino-like, respectively. 

Finally, we consider two cases of heavy (600 GeV) 
neutralinos. The first (gaugino-like) annihilates mostly 
(87%) to bb and to t + t~ or tt the remaining times. 
The second (higgsino-like) annihilates 65% of the time 
to gauge bosons, 19% to t+t~ and 14% to higgs bosons. 

Although these five cases do not completely encompass 
the very wide array of characteristics neutralinos may 
have, they do describe effective benchmarks. Further- 
more, the relevant characteristics of a LSP with a mix- 
ture of these annihilation modes or another mass could be 
readily inferred by inspecting these representative mod- 
els. 

The effect of propagation on the positron spectrum 
depends strongly on the distance from the source (the 



dark matter clump) and the observer. Since the tidal 
radius of typical clumps is typically much smaller than 
the distances to clumps we will consider, we can treat 
such a clump spatially as a point source. 

To compare to the data recorded by HEAT, rather 
than simply the positron flux, we consider the ratio of 
the positron flux to the combined positron and electron 
fluxes, called the positron fraction. We do this by using 
the spectra for secondary positron, secondary electrons 
and primary electrons found in Ref. [33] . 

Figures 3-7 show the positron fraction, as a function 
of positron energy, for each of the five representative su- 
persymmetry scenarios described above. We varied the 
distance to the source (the dark matter clump) to deter- 
mine the effect on the positron spectrum. In each figure, 
the solid line represents the distance at which the pre- 
dicted spectrum is best fit to the data. In all five cases, 
we found very good fits to the HEAT observations. 

Dotted and dashed lines represent the spectra for a 
source at (1 and 2a) distances less than and greater than 
found for the best fit, respectively. For these lines, the \ 2 
is larger by 1 and 4, respectively (la and 2a). The nor- 
malization was considered to be a free parameter. The 
predicted spectrum is compared to the error bars of the 
94-95 and 2000 HEAT data. These results are summa- 
rized in tables I and II. 

Note that our results are quite different from that of 
other collaborations which treat the spatial distribution 
of WIMP annihilations as constant. Introducing a spatial 
distribution of clumps allows the data to be well fit by our 
predictions for SUSY models which would not otherwise 
be consistent. 
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FIG. 3. The positron fraction, as a function of positron 
energy (in GeV), for a 50 GeV neutralino which annihilates 
96% to bb and 4% to t + t~. The solid line represents the 
distance to the dark matter clump at which the predicted 
spectrum best fits the data. Dotted and dashed lines represent 
the spectra for a source at (1 and 2a) distances less than 
and greater than found for the best fit, respectively. The 
normalization was considered to be a free parameter. The 
error bars shown are for the 94-95 and 2000 HEAT flights. 
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FIG. 4. The predicted positron fraction, as a function of 
positron energy (in GeV), for a 150 GeV neutralino which 
annihilates 96% to bb and 4% to r + r _ . Otherwise, the same 
as in figure 3. 
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FIG. 6. The predicted positron fraction, as a function of 
positron energy (in GeV), for a 600 GeV neutralino which 
annihilates 87% to bb and to t + t~ or it the remaining times. 
Otherwise, the same as in figure 3. 
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FIG. 5. The predicted positron fraction, as a function of 
positron energy (in GeV), for a 150 GeV neutralino which 
annihilates 58% to W + W~ and 42% to ZZ. Otherwise, the 
same as in figure 3. 
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FIG. 7. The predicted positron fraction, as a function of 
positron energy (in GeV), for a 600 GeV neutralino which 
annihilates 65% to gauge bosons, 19% t + t~ and 14% to Higgs 
bosons. Otherwise, the same as in figure 3. 
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TABLE I. The distances to a dark matter clump found to 
best fit the observed positron spectrum for the five supersym- 
metry models described in the text. One and two sigma error 
bars are also given. Cases 1, 2 and 4 correspond to models 
dominated by annihilation to bb (with a small fraction anni- 
hilating to t~t + or, in the heaviest case, tt). Cases 3 and 5 
correspond to models which annihilate dominantly into gauge 
bosons (with, in the heavy small fraction annihilating 

to bb, t~t + or it). 



TABLE II. The rate of neutralino annihilations (per sec- 
ond) in a dark matter clump required to account for the ob- 
served positron excess, using the distances in the previous ta- 
ble. Again, the five supersymmetry models which have been 
described in the text are considered. Results are shown for a 
clump the best-fit distance from Earth, and for a clump at a 
distance corresponding to 2a from the best fit. 
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VI. ASSESSMENT 

The results shown in figures 3-7 and tables I and II can 
be used to help assess the size and density profile of a 
dark matter clump required to account for the observed 
positron excess. The annihilation rate in a clump is a 
function of the annihilation cross section, the LSP mass 
and clump characteristics (mass and density profile): 

R =\ <°rf™ ,,, (4, 

LSP ^clump 

where < av > is the low-velocity neutralino annihilation 
cross section times velocity, tolsp is the neutralino mass, 
■^ciump is the mass of the clump, V c i ump is the volume 
of the clump, and / is the enhancement from uneven 
distributions of mass within the clump, defined in section 
IV. 

The first factor in the expression for the annihila- 
tion rate depends on particle physics. To ascertain how 
large this factor (< av > /ni 2 SP ) could be, we ran 
a Monte Carlo, randomly selecting parameters of the 
MSSM. This method found several thousand SUSY mod- 
els which respected all collider constraints and provided 
a neutralino LSP with a relic density in the range mea- 
sured by WMAP (0.095 < O^o/i 2 < 0.129). The models 
were selected by randomly varying values of Mi between 
10 and 10,000 GeV, |/z| between 10 and 10,000 GeV (with 
either sign), the mass of the pseudoscalar Higgs, m^, be- 
tween 30 and 1000 GeV, tan (3 between 1 and 60, and the 
squark and slepton masses between 30 and 10,000 GeV. 
The trilinear couplings, A t and A^, were varied between 
30 and 10,000 GeV, as well. M Y and M 3 were deter- 
mined by the GUT relations. Of the models our Monte 
Carlo produced, none with LSP masses in the range of 
30-70 GeV had values of < av > /ttjlsp larger than 
about 2 x 10~ 29 cm 3 s _1 GeV -2 . For the other four char- 
acteristic SUSY scenarios we are considering through- 
out this paper, we find upper values for this factor of 
4 x lO- 30 (150 GeV, gaugino), 2 x lO" 30 (150 GeV, hig- 
gsino), 1 x 10~ 31 (600 GeV, gaugino) and 1 x 10~ 31 (600 
GeV, higgsino) cm 3 s _1 GcV~ 2 . We will use these values 
for < av > /m-Lgp for the remaincr of our calculation in 
an effort to be optimistic. 

For a light LSP (m x o < m w ±), annihilations will be 
dominated by bb except for a small fraction of models 
which have very small annihilation cross section and are, 
therefore, not of interest to us. As we said before, for 
a 50 GeV neutralino, < av > /m 2 SP can be as large 
as - 2 x 10~ 29 cm 3 s" 1 GeV -2 . Combining this with the 
result of table II, this leads to the requirement of a clump 
with / x M 2 ump /U clump ~ 2.3 x 10 7 M|/pc 3 using the 
best fit distance and 2.0 - 5.0 x 10 7 M 2 /pc 3 using the 
distances corresponding to a 2a fit. These results and 
the results for our other supersymmetry benchmarks are 
shown in table III. 



Now that we have determined the nature of the clump 
required to account for the observed positron excess, we 
can begin to estimate the probability of such a clump 
being present. For example, considering a 50 GeV neu- 
tralino, if we require a clump within the best-fit distance 
of 0.27 kpc and with at least / x M^ lu /V^i ump ~ 2.3 x 
10 7 M,|/pc 3 , using the results of figure 2, we see that the 
probability of this occuring is approximately the number 
density of such clumps multiplied by the volume within 
this distance, (2 - 3 x 10" 4 ) (4tt x 0.27 3 /3) ~ 2.5 x 10~ 5 , 
or one in 40,000. These results and the results for other 
SUSY models are shown in table IV. 

Note that if non-thermal processes are responsible for 
the observed density of dark matter [34], the neutralino 
annihilation cross section may be considerably larger and 
our estimates inaccurate. 
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TABLE III. The values of / x M 2 lump /V c i ump , in units of 
Mq /pc 3 , required from a dark matter clump to account for 
the observed positron excess. The five supersymmetry mod- 
els described in the text are again used. These values are 
given for the distance to the clump which results in the best 
fit spectrum as well as for the distances corresponding to 
the two sigma error bars around the best fit. The values of 
< av > /ttilsp used are fixed to optimistic values of 2 x 10 -29 , 
4 x 10~ 30 , 2 x 10~ 30 , 1 x 10~ 31 and 1 x 10~ 31 cm 3 s" 1 GeV -2 , 
for each of the five supersymmetry models used, respectively. 
Cases 1, 2 and 4 correspond to models dominated by annihi- 
lation to bb (with a small fraction annihilating to t~t + or, 
in the heaviest case, tt). Cases 3 and 5 correspond to models 
which annihilate dominately into gauge bosons (with, in the 
heavy small fraction annihilating to bb, r r + or tt). 



SUSY Model 


Usinj 


; Best-Fit Dist. 


Using 2a Dist. 


m x o=50 GeV (Case 


1) 


3 x 10~ b 


5 x 10~ 4 


m x o=150 GeV (Case 


2) 


1 x 10~ 4 


8 x 10~ 4 


m x o=150 GeV (Case 


3) 


2 x 10~ b 


2 x 10~ 4 


m x o=600 GeV (Case 


4) 


£ io~ b 


- 10~ 4 - 10~ b 


m x o=600 GeV (Case 


5) 


< 10" b 


~ 10" 4 - 10" b 



TABLE IV. The probability of a clump with sufficient 
/ x M 2 lump /V c iump within the best fit and 2a distances. Other 
details are the same as in table III. The most optimistic val- 
ues for cross sections and substructure number densities were 
used. 
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VII. CONCLUSIONS 

In each characteristic supersymmetry model consid- 
ered, we found that (up to normalization) the spectrum 
of positrons observed by HEAT could be well fit by a 
dark matter clump at an unknown distance acting as a 
point source of positrons. Furthermore, by treating dark 
matter clumps as point sources rather than using an even 
distribution, we can estimate the distance to the source 
clump for a given neutralino mass and dominating anni- 
hilation modes. Such distances are typically 0.2 to 2.0 
kiloparsecs, considerably larger than the average size of 
the core or tidal radius of such clumps. Thus the point- 
like approximation for annihilations is justified. In the 
unlikely case that the solar system resides inside of such 
a clump, the spectrum will be quite sharply rising at high 
energies, fitting the data poorly, as can be inferred from 
figures 3-7. 

Once a source distance is estimated, we can determine 
the annihilation rate which must occur in the clump 
to properly normalize the positron flux. This rate is 
the result of both the characteristics of the neutralino 
(< (tv > /m£ SP ) and the clump (/ x M^ ump /V c i ump ). 
By considering a model of galactic substructure which 
can predict the number density of dark matter clumps, 
we can use the required annihilation rate to estimate the 
probability of a clump being present which is sufficient 
to produce the observed positron excess for a given su- 
persymmetry model. It is likely that our substructure 
model overestimates the number of clumps. With this in 
mind, we consider this as a rough upper limit. 

We find that, using optimistic neutralino annihilation 
cross sections and annihilation modes, the probability of 
local dark matter substructure being sufficient to provide 
the observed positron excess is quite small, on the order 
of 10~ 4 or less. Given this, we conclude that although it 
would not be impossible for annihilating supcrsymmctric 
dark matter to produce the observed positron excess, it 
is a very unlikely and unnatural solution. 

Acknowledgments: DH and JT are supported by the 
Leverhulme Trust. 
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